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ABSTRACT 

We investigate the large-scale angular distribution of the short-Gamma ray bursts (SGRBs) 
from BATSE experiment, using a new coordinates-free method. The analyses performed take into 
account the angular correlations induced by the non-uniform sky exposure during the experiment, 
and the uncertainty in the measured angular coordinates. Comparising the large-scale angular 
correlations from the data with those expected from simulations using the exposure function 
we find similar features. Additionally, confronting the large- angle correlations computed from 
the data with those obtained from simulated maps produced under the assumption of statistical 
isotropy we found that they are incompatible at 95% confidence level. However, such differences 
are restricted to the angular scales 36° — 45°, which are likely to be due to the non- uniform 
sky exposure. This result strongly suggests that the set of SGRBs from BATSE are intrinsically 
isotropic. Moreover, we also investigated a possible large-angle correlation of these data with the 
supergalactic plane. No evidence for such large-scale anisotropy was found. 

Subject headings: large-scale structure of universe — gamma rays: bursts — methods: statistical 



1. Introduction 

The apparent isotropy in the large-scale an- 
gular distribution of the Gamm a ray bursts 
(GRBs) is a long-standing debate (iMeeean et al 



19981: iMeszaros. Bagolv. fc Vavrekl l2000t )._Since 



1992t iBriggs et al.l [l99a iTegmark et all Il996t 



Piran &: Singhll997trMetzger et aljll997t Balazs. Meszaros^^^ l lorva'^h 



by cross-correlating them with catalogs of cos- 
mic objects, e. g., early- type galaxies, hard X-ray 
sources, etc. ('Briggs et al.' 'l996'; iTegmark et al.l 
11996; Piran 2004; M cszar os 2006). 

The reported statistical analyses of the all- 
data from BATSE show that their 
liffUfe-scai e' angular distribution is cons istent with 



,1 n ; , , 7- rrr^ ,r^r A „•/ ^77~^ Tl — '. — T^otropy dPiranl [200I; iMeszaroj l2006l ). although 

the farst detection with the VELA satellite ( Klebesadel et al.1 ^ ,. 1 ,•<■ ^ , , , 

, . . p , , . 1 1 . ^ aspects like observational artifacts have not been 



19731 ) the origin of these highly energetic events 
has remained a challenge. Even if the origin 
of GRBs turns out to be extra-galactic or cos- 
molo gical, as sugp :es ted by current da t a (see , 
e.g.. IPiranI ('20041: ' Zhang fc Meszarod (|2006t ): 
MeszarosI (|2002. ,200^ and references therein) , 
this does not ensure that their distribution is 
isotropic. Up to now, no dominant anisotropics 
has been found in the angular distribution of 
GRBs. However, if detected, small anisotropic 
effects may reveal valuable information about 
their origin. Additionally, the discovery of a large 
angular scale pattern in the sky distribution of 
GRBs may be useful to identify their sources 



fully considered in some of these studies. It is well 
known that anisotropics with distinct origins man- 
ifest themselves on different angular scales and 
with different magnitudes. In this connection, it 
is reasonable to consider different approaches that 
can, in principle, provide information about multi- 
ple types of anisotropy, imprinted as angular cor- 
relation signatures (ACS), that may be possibly 
present in the data. 

Here we apply a new coordinates-free method 
to search for large-scale AC S (> 18°) in a subse t 
of the BATSE GRB data (jMeegan et al.ll2000l ). 
namely the Short-GRBs, and then investigate 



their significance levels through the comparison 
with a large number of Monte Carlo maps. Such 
simulated maps were produced under similar con- 
ditions as the catalog under analysis, that is, tak- 
ing into account the non-uniform sky exposure of 
BATSE and the uncertainty in the coordinates 
measurements. Furthermore, for completeness, 
we also compare the ACS of the catalog of GRBs 
with those corresponding to statistically isotropic 
Monte Carlo maps. Finally, we also investigated 
the possible large-scale angular correlation be- 
tween the set of Short-GRBs and the supergalac- 
tic plane, in an attempt to search for likely host 
ga laxies of these events (as recently suggested 
bv lGhirlanda et all (l2006h ). 

The outline of this paper is the following: in 
section [2] we use GRBs data from the BATSE ex- 
periment to determine the Short-GRBs catalog to 
be investigated, and in section [3l we describe the 
method employed in such studies. The data anal- 
yses and results are shown in section |4l and finally 
in section [5] we formulate our conclusions. 

2. The Short-GRBs from BATSE catalog 

The physical analysis of GR Bs utilizes their 
temp oral and spectral properties (jFishman fc Meegan 



19951 ). Despite the different light-curves observed 
in the spectra of GRBs, a useful parameter to 
classify them is the burst duration Tgo, defined as 
the time interval during which 90% of the fluence 
is measured. The current BATSE catalog 4Bc 
contains 2 702 events, where onl y 2 037 GRBs 



have their parameter Tgo measured (jMeegan et al 
l200fl[ ). 

At first, the Tgo value was used to divide the set 
of GRBs into two different sub-classes: the Short- 
GRBs (SGRBs), with Top < 2s, and the Long 
GRBs, with Tgo > 10s iKouveliotou et al. 199 



Zhang fc Meszdrod l2006t iMeszarosI I2002L 1200 



However, the use of this definition of SGRBs is 
instrument depen dent and is susceptib le to ob- 
servational biases ( Hakkila et al. 2007b[) . For this 
reason, one should consider in a ddition to the Tgp 



criterium the parameter if i?3/2i (iMukherjee et al 



1998D which is defined as the 100 to 300 keV 



fluence divided by the 2 5 to 100 keV fluence of 
each GRB in BATSE 4Bc (|Hakkila et al.ll2007allbl ). 



Thus, the a ppropriate defi nition for a catalog of 
SGRBs is fflakkila et all l2007allbh : C = {516 
events with 2s < Tgo < 4.7s and HR3/21 > 3}. 

With this information, and using a new coordinates- 
free method to be described in the next section, 
we shall perform a detailed analysis of the large- 
scale ACS present in the sky distribution of the 
SGRBs from BATSE. 

3. The 2PACF and the Sigma-Map analy- 
sis 

Let f2j° = Q{0j,(j)j;^o) G <S^ be a spheri- 
cal cap region on the celestial sphere, of 70 de- 
grees of aperture, with vertex at the j-th pixel. 



] 



1, 



, TVcaps, where 



4>j) are the angular 
coordinates of the center of the j-th pixel. Both, 
the number of spherical caps A^'caps and the coor- 
dinates of their center {dj,4>j) are defined using 



the H EALPix pixelization scheme (|G6rski et al 
20051 ). The spherical caps are such that their 



union completely covers the celestial sphere S^. 

Let be the catalog of cosmic objects located 
in the j-t h spherical cap f^J ° The 2PACF of these 
objec ts ( Chen fc Hakkilal Il998l : iPadmanabhan 
19931 ). denoted as Tj(7i;7o), is the difference 
between the normalized frequency distribution 
and that expected from the number of pairs- 
of-objects with angular distances in the inter- 
val (7i - 0.55, 7j -I- 0.5(5], i = 1, . • . , A^bins, where 
7i = (i — 0.5)6 and S = 270/A^bms is the bin- width. 
The expected frequency distribution is achieved 
by a large number of Monte Carlo realizations of 
isotropically distributed objects in contain- 



ing a si milar number of obje cts as in C-' (jTeixeira 
2003; B ernui fc Villelall2006l ). The 2PACF has the 
property that its mean, obtained by integratin g 
over all separation angles (jChen fc Hakkila 1998[ ). 



is zero. A positive (negative) value of Tj indi- 
cates that objects with these angular separations 
are correlated (anti-correlated), while zero indi- 
cates no correlation. 

Define now the scalar function ct : fij 1-^ 3?^, 
for j = 1 , . . . , Acaps J which assigns to the j- 
cap, centered at {6j,(j)j), a real positive number 
ffj = a{9j,4>j) e 5R^. The most natural way of 
defining a measure a is through the variance of the 



^ http: / / gammaray. msfc. nasa. gov /batse /grb/ catalog/current / 



2 



Tj function, we thus define (jBernui et al.l 120071 ) 



^bins 

E 



Tj(7i;7o)- 



(1) 



To obtain a quantitative measure of the ACS of 
the GRBs sky map, we cover the celestial sphere 
with A^caps spherical caps, and calculate the set of 
values {aj, j = 1, TVcaps} using eq. (P). Patch- 
ing together the set {aj} in the celestial sphere 
according to a coloured scale (where, for instance. 



, minimum 



blue, cr™'^«mum ^ ^g^) qI^^j^Jj^ 

a sigma-map. Finally, we quantify the ACS of a 
given sigma-map by calculating its angular power 
spectrum. Because the sigma-map assigns a real 
value to each pixel in the celestial sphere, that 
is (J = a{9,(j)), one can expand it in spherical 
harmonics: (t(6', (^) = J^e.m ^em^e^i^' 'f') where 
the set of values {Si, £ = 0,1,2,...}, defined by 



Si ^ (1/(2£+1))E: 



em\ 



is the angular 



power spectrum of the sigma-map. Because we 
are interested in the large-scale angular correla- 
tions, we shall concentrate on {Se, i — 1,2, 10}. 
Notice that we are interested in the angular power 
spectrum of the sigma-map, and not that of the ce- 
lestial sphere where the GRB events are located; 
this la ter c ase was already done by iBriggs et al 



(|l996l ) and iTegmark et"aLl (|l996f) . As we shall 
see, the sigma-map analysis is able to reveal very 
small anisotropics, like those induced by the BAT- 
SEs sky exposure, despite the small burst detec- 
tion rate of the SGRBs. 

4. Data analyses and results 

In this section we shall apply the sigma-map 
method explained in the previous section to study 
the large-scale ACS present in the angular distri- 
bution of the 516 SGRBs listed in the catalog C. In 
the following, all the sigma-maps were calculated 
using spherical caps of 70 = 90° of aperture, that 
is hemispheres (smaller spherical caps have less 
SGRBs in each , j — I, ■ ■ ■ , Acaps, hence pro- 
duce large statistical noise in the Tj functions). 
We also used A^bins — 90 and A^caps = 768 in these 
analyses. 

An important issue that deserves close inspec- 
tion is the presence of anisotropic ACS in the data 
induced by the non-uniform sky exp osure (NUSE) 
during the BATSE experiment (jHakkila et al 



Il998l ). expected because some latitudes of the 
sky were over-observed while others were under- 
observed. Because there are no reports quantify- 
ing or tracing out the infiuence of the NUSE at 
large ang ular scales in the c urren t BATSE catalog 
4Bc (see IChen fc Hakkilal ()l998f ) for analyses of 
the 3B and 4B catalogs) it is interesting to use 
the sigma-map method to investigate the possi- 
ble anisotropic angular correlations that may be 
present in the data even if their magnitudes are 
small. For this, our strategy to reveal the large- 
scale ACS in the data runs in three steps. First, 
we produce 1 000 Monte Carlo maps simulating 
the sky positions of 516 co smic objects accord - 



ing to the NUSE function (jHakkila et al.lll998l ). 



then we calculate in each case the corresponding 
sigma-map, and finally we compute the angular 
power spectrum {Si, £ — 1, 10} of each of these 
sigma-maps. Second, we generate 1 000 Monte 
Carlo maps simulating the sky positions of 516 
isotropically distributed cosmic objects, then we 
compute for each of these Monte Carlo maps their 
corresponding sigma-maps, and finally we calcu- 
late the angular power spectrum of these sigma- 
maps. Third, we calculate the sigma-maps, and 
their respective angular power spectrum, of the 
SGRBs listed in catalog C. 

In figure [1] we show two sigma-maps in galac- 
tic coordinates. In the top panel, we show the 
average of 100 sigma-maps, randomly chosen in 
between the 1 000 sigma-maps computed from a 
similar number of Monte Carlo sky maps which 
simulate different catalogs of SGRB according to 
the NUSE function. In the bottom panel we ex- 
hibit the sigma-map corresponding to the catalog 
C. 

In figure [2] we display a comparative analyses, 
taking into account isotropic and non-isotropic 
cases, of the angular power spectrum of the sigma- 
map obtained from the angular distribution of the 
SGRBs listed in C. In the top panel, we plot- 
ted the angular power spectrum Se versus £ of the 
sigma-map obtained from the catalog C, together 
with the mean of 1 000 sigma-maps computed from 
an equal number of statistically isotropic Monte 
Carlo sky maps. In the bottom panel, the plot is 
similar except that the Monte Carlo sky maps have 
anisotropic ACS because were produced consider- 
ing the NUSE function of BATSE. In both plots 
the shadowed areas correspond to 2 standard de- 
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viations. Besides some small differences, the an- 
gular power spectra corresponding to the sigma- 
maps computed from the SGRBs show a very sim- 
ilar large-scale structure when compared with the 
mean angular power spectrum of the sigma-maps 
obtained from Monte Carlos produced according 
to the NUSE function. 

A comparative analysis of the ACS correspond- 
ing to these cases, isotropic and non-isotropic 
due to NUSE function, is better seen if we plot 
i (^+1) Si versus ^, as showed in figure[3] There we 
display the correponding data from the SGRBs to- 
gether with the mean of the angular power spectra 
of the isotropic and non-isotropic cases, where now 
the shadowed area corresponds to 2 standard devi- 
ations of the isotropic case. As observed, the data 
have a very similar behavior to the non-isotropic 
case, and is different from the isotropic case which 
shows a flat spectrum. Thus, data and simulated 
isotropic maps are incompatible at 95% confidence 
level. However such differences are mainly re- 
stricted to the angular scales 36° — 45° which ex- 
actly reproduce the imprints exhibited by the an- 
gular power spectrum of the non-isotropic case. In 
the absence of ACS other than those expected by 
the NUSE of the BATSE experiment, this resuh 
strongly suggests that the SGRBs are intrinsically 
isotropic. 

To test the robustness of our calculations 
we also performed the sigma-map analyses with 
A^bins = 180 and Acaps = 3 072, obtaining the 
same result. 

Furthermore, we also searched for a possible 
correlation between the SGRBs listed in C with 
the supergalactic plane, where nearby galaxies ap- 
pear to be more concentrated. Because we do 
not know how many events could be originated 
in these galaxies, we generate three sets of 300 
Monte Carlos considering in each case a differ- 
ent number of simulated GRBs provided by an 
anisotropic distribution which selects events near 
the supergalactic plane. That is, we generated 
sets of maps where 33%, 50%, and 66% of the 
events were produced by such anisotropic distri- 
bution, respectively. We then computed their cor- 
responding sigma-maps in order to compare them 
with the sigma-map calculated from the SGRBs 
listed in C. To measure such a possible correlation 
we computed the linear Pearson correlation coef- 
ficient between the sigma-map of the SGRBs and 



each one of the sigma-maps obtained from these 
sets of Monte Carlo realizations. Notice that such 
a coefficient varies from (for totally uncorrelated 
maps) to 1 (for fully correlated maps). Our results 
show that the Pearson's coefficient is, in mean, less 
than 0.03 (using 300 Monte Carlos for each of the 
three above mentioned cases). To realize whether 
this value is statistically significant, we computed 
the Pearson's coefficient in some illustrative cases. 
For example, the mean Pearson's coefficient corre- 
lating one sigma-map, coming from a given set of 
sigma-maps computed from the above mentioned 
66% anisotropic Monte Carlo maps, with the rest 
of sigma-maps from such set is 0.23 ± 0.16 (the 
result is similar in the other two cases). On the 
other hand, the mean Pearson's coefficient corre- 
lating one sigma-map, chosen randomly from the 
set of 1 000 sigma-maps calculated from Monte 
Carlo maps produced according to the NUSE func- 
tion, with the rest of sigma-maps of this set is 
0.25 ± 0.16. Similarly, the mean Pearson's coef- 
ficient correlating any sigma-map, from the set 
of 1 000 sigma-maps computed from Monte Carlo 
statistically isotropic maps, with the rest of sigma- 
maps of this set is 0.15 ± 0.12. 

Moreover, comparing the sigma-map computed 
from BATSE SGRBs catalog C with each of the 
1 000 sigma-maps, obtained from Monte Carlo 
maps generated according to the NUSE func- 
tion, the mean Pearson's coefficient is 0.14 ±0.11. 
A similar analysis of the sigma-map of BATSE 
SGRBs but now considering those Monte Carlo 
maps generated under the statistical isotropy hy- 
pothesis, tells us that the Pearson's coefficient is, 
in mean, 0.12 ±0.09. Taken together this informa- 
tion suggests that there is no evidence for a large- 
scale angular correlation between BATSE SGRBs 
and simulated maps produced considering differ- 
ent amounts of events coming from an anisotropic 
distribution that selects positions in the super- 
galactic plane and its surrounds. Notice that this 
res ult does not contradic t the correlation found 
bv lChirlanda et aL ( 2006f ) which is valid for small 
angular distances < 3°, while the present analysis 
concerns angular scales > 18°. 

Finally, we also tested the robustness of our re- 
sults under the change of the angular coordinates 
of the BA TSE SGRBs due to the measured er- 
ror boxes ( Briggs et al. 19981 ). This was done by 



sorting their angular coordinates within the lim- 
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its given by the error boxes (±2.5°) we found no 
measurable difference with respect to the results 
presented here. 

5. Conclusions 

The purpose of this study is to know the large- 
scale angular correlations of the set of 516 BATSE 
SGRBs, and to discover if these correlations are 
compatible with a statistically isotropic distribu- 
tion of events, or instead they reveal the ACS re- 
sulting as a consequence of the NUSE function of 
BATSE experiment. To elucidate this, we need 
to know the angular power spectra of two sets of 
sigma-maps: one set is computed from statistically 
isotropic Monte Carlo maps and the other is cal- 
culated from Monte Carlo maps that simulate the 
sky position of the events using the NUSE function 
of BATSE experiment. After that, we compare the 
power spectra of these sigma-maps with the angu- 
lar power spectrum of the sigma-map computed 
from the BATSE SGRBs data. 

The first thing to be noticed in the angular 
power spectrum of the BATSE SGRBs, plotted in 
figure[2l is the absence of dominant anisotropics at 
the largest angular scales 60° - 180° {£ = 1,2,3), 
and a similar situation for angular scales < 30° 

> 6). However, peculiar features appear at the 
angular scales 36° — 45° {£ = 4, 5), and for a bet- 
ter understanding of what information is encoded 
there we plot these data in the form £{£ + 1) Si 
versus i (see figure [31). In fact, figure [3] reveals a 
second interesting thing, that is, the non-flat spec- 
trum of the BATSE data (represented by bullets) 
which clearly differs from the flat angular power 
spectra showed by the statistically isotropic Monte 
Carlo data (the dashed line). We also observe in 
figure [3] that the mean of the angular power spec- 
tra of the sigma-maps computed from Monte Carlo 
maps generated according to the NUSE function 
(the dot-dashed line) has also a non-flat spectrum 
which is very similar to the corresponding one 
obtained from BATSE SGRBs. In other words, 
the large-scale angular correlations of the BATSE 
SGRBs exhibit the anisotropic imprints expected 
in the data due to the NUSE of BATSE experi- 
ment. Other ACS are not found to be statistically 
significant, at 95% CL. In conclusion, these results 
strongly suggests that the SGRBs are intrinsically 
isotropic. 



Finally, we also studied the possible large- 
angle correlation between the SGRBs data and 
Monte Carlos with (different amounts of) sim- 
ulated events concentrated towards the super- 
galactic plane. No evidence for such large-scale 
anisotropy was found in the BATSE SGRBs. 
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Bottom: This is the sigma-map calculated from 
the catalog C of BATSE SGRBs data. 
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Fig. 2. — Averaged angular power spectra of 
1 000 sigma-maps computed from equal number 
of Monte Carlo maps, which were produced con- 
sidering two cases: isotropic (top) and anisotropic 
(bottom) hypotheses. For the anisotropic case, the 
Monte Carlo maps were generated using the NUSE 
function. Together we plot the angular power 
spectrum of the sigma-map, represented as bul- 
lets, corresponding to the BATSE SGRBs listed in 
the catalog C. The shadowed regions correspond 
to 2 standard deviations in each case. 
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Fig. 3. — Plot of the angular power spectra in the 
form £ (^+1) versus £. The dashed (dot-dashed) 
line corresponds to the mean angular power spec- 
trum of 1 000 sigma-maps, each one computed 
from a Monte Carlo isotropic map (anisotropic 
map, according to the NUSE function). The shad- 
owed region corresponds to 2 standard deviations 
of the isotropic case. Together with these data we 
plotted the angular power spectrum of the SGRBs 
listed in C, represented as bullets. 



